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Abstract 

A theoretical model of parametric magnetostrictive generator of ultrasound is considered, taking into account magnetic 
and magnetoacoustic nonlinearities. The stability and temporal dynamics of the system is analized with standard techniques 
revealing that, for a given set of parameters, the model presents a homoclinic or saddle-loop bifurcation, which predicts that 
the ultrasound is emitted in the form of pulses or spikes with arbitralily low frequency. 

PACS numbers: 05.45.-a, 75.80.+q, 43.25.+y 



I. INTRODUCTION 

Magnetostriction is essentially a nonlinear phe- 
nomenon, which accounts for the change of dimension of 
a magnetic material under applied magnetic fields. How- 
ever, most of the applications of this phenomenon con- 
sider small amplitude fields and consequently exploit only 
their linear properties, as is the case of electromagnetic- 
acoustic transducers. The consideration of nonlinearities 
in the description of this process lead to the discovery of 
new properties, such as the parametric wave phase conju- 
gation (WPC), a phenomena which is actually an active 
field of research with promising applications into acous- 
tic microscopy p] or harmonic imaging 0. A review on 
WPC theory and methods can be found in 0. There 
is also increasing interest in the developement of mag- 
netostrictive transducers working at high powers, where 
nonlinear effects are not negligible. The advances in this 
field come in paralel with the search of novel magnetic 
materials with high magnetostriction values. 

On the other side, parametric phenomena in different 
fields of nonlinear science share many common features, 
such as bistability, self-pulsations and chaos among oth- 
ers. Is precisely this analogy that motivates our search for 
complex dynamical phenomena in parametrically driven 
magnetoacoustic systems. 

Different models have been proposed for the descrip- 
tion of magnetoacoustic interaction in ferromagnetic ma- 
terials. Also, many experimental progress in this field has 
been achieved, and the main parameters involved in the 
process have been determined 4] . In this paper we anal- 
ize the dynamical behaviour of this system in the pres- 
ence of magnetic nonlinearity. It is shown that a cubic 
magnetic nonlinearity is responsible for the appearance 
of new effects not reported before in this system, such as 
sclf-pulsing dynamics and spiking behaviour, related to 
the existence of homoclinic bifurcations 



II. THE MODEL 

The physical system considered in this paper consists 
in a electric RLC circuit, driven by an external ac source 
at frequency 2u and variable amplitude e. The induc- 



tance coil, with density of turns n, transverse section S 
and length L, contains a ferromagnetic ceramic material 
which acts as an acoustical resonator and is the origin of 
nonlinearities in the system when the driving e is high 
enough. A theoretical model for this system has been 
derived in the resonant case in We review here the 
details of the derivation for convenience. 

The equation for the circuit is er + ec + £z = text, 
where e ext = £cos(2wi), er — IR, ec — q/C, and el 
depends nonlinearly on the magnetic fields. 

We consider two effects that introduce the nonlinear- 
ities in the problem. The first source of nonlinearity re- 
sults from magnetoacoustic interaction, which appears 
via phonon-magnon proccesses. In the case of parametric 
magnetostriction ceramics, the generation of ultrasound 
takes place at half of the frequency of the driving. This 
process is represented, in general, by a lagrangian den- 
sity term in the form Li nt = aijkHiUjUk,\H$ where Uk is 
the acoustic displacement component in the direction k 
(k = x, y, z) and a has in general tensor character. 

The total magnetic induction H acting along the axis 
of the material is the result of three contributions: a 
static field Hq produced, e.g. by a permanent magnet 
surrounding the active ceramic material or a coil carrying 
a stationary current, an alternating field H q (t) induced 
by the ac current in the circuit, and finally a field pro- 
duced by material deformations Hi nt , which results from 
the magnetoacoustic interaction. From the lagrangian 
density it follows that, in the case of longitudinal waves 
propagating in the z axis (u x = u y = 0), the acousto- 
induccd magnetic field has the form 




where a = a zzz is the coefficient of coupling between the 
active medium of volume V and the pump source. Thus, 
the effective magnetic induction takes the form 

H = H + H q (t)+H int (t). (2) 

where we assume that the relation Hq >> H q (t), Hi nt (t) 
holds. Alternatively, H given by ^ can be considered 
as a nonlinear pumping field 6]. Note that the appear- 
ance of the subharmonic field saturates the magnetic field 
value, and modifies the effective pumping. 



A second source of nonlinearity is typical of ferromag- 
netic materials. We assume that, for weak saturation, the 
scalar nonlinear relation between the magnetic field and 
the magnetic induction, B = /i (H) H, can be written, 
following 7], as 



(3) 



where /i is linear the permeability of the material and 
the third order magnetic susceptibility, which in turn 
depends on the frequency. 

Applying the Faraday law under the previous assump- 
tions, considering only resonant terms oscillating at the 
frequency of the driving 2lo, and neglecting terms higher 
than quadratic in H q and Hi nt , we get 



(4) 



d f d f 

UN— / H mt ds + fi oX {3) N— / H H q H int ds, 

where / = dq/dt is the electrical current in the circuit, C 
is the inductance and H q = NI/L. In terms of the charge 
in the capacitor, and taking into account that V = SL 
and that n — N/L is the density of turns, the circuit 
equation takes the form 
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where the last two terms represent the nonlinearities re- 
lated magnetoelastic interaction and magnetic nonlinear- 
ity respectively. 

Let us consider now the evolution of the acoustic wave. 
Considering the current induced magnetic field H q as the 
main source term, we find |6| 



o Trpr - V 2 u = aHJt)u = an^-u, 
v z ot z 



dq 
dt 
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where v is the propagation velocity of sound in the mate- 
rial. Here the effect of the acousto-induced magnetic field 
(quadratic in the small coupling constant a) has been ne- 
glected. However, we have checked that the consideration 
of this small term do not imply qualitative changes in the 
results reported below. 

We consider solutions of Eqs. JSJ and © in the form of 
quasi-harmonic waves, i.e., whose amplitudes are slowly 
varying in time. In this case we can write 



q (t) = i [Q (t) exp (2iwi) + c.c] , 
u (r, t) = i [U (t) exp (iujt) + c.c] g (r±) sin (kz) 



(7a) 

I, 
(7b) 



where u (r, t) has the form of cavity modes of the acous- 
tical resonator, k = rmr/L define the cavity resonances, 



2u> = 1/vxC is the pumping frequency resonant with 
the circuit, and 



dX 



dt 



« \loX\ 



(8) 



with X = Q or U, represents the slowly- varying envelope 
approximation. 

Under these assumptions, the slow amplitudes obey 
the evolution equations 



dQ 

dt 
dU 

~dt 



= tE- lQ Q + Z 1 U 2 + i£ 2 \U\ 2 Q 1 (9a) 
= -1uU + UQU*> (9b) 



where E = £ /AojC, and 7q = R/2C and 7^ represent the 
electric and acoustic losses respectively. The last param- 
eter is introduced phcnomcnologically, and take into ac- 
count the losses due mainly to radiation from the bound- 
aries. The coefficients of the nonlinear terms are defined 
as 



£1 = 
£3 = 



UnaL f 2 
-^]9{^) dr„ 



(10) 



Equations 10 are the model considered in |5| with some 
corrections. This system can be further simplified using 
the normalizations 



which leads to the model 

X + Y 2 + i V \Y\ 2 X, 
XY*), 



dX 



(12) 



dY 



where 1 — lu llQ ' IS the ratio between losses, r = jgt is 
a dimensionless time and 



V 
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(13) 



remains as the single parameter of nonlinearity, with 
fi r = /i//i . Due to the number of parameters involved 
in 77, it can be varied over a wide range of values. Note 
that Eqs. Ijl2(l also possess the symmetry Y — ► — Y. 



III. STATIONARY SOLUTIONS AND STABIL- 
ITY 

Equations (|12fl possess two kinds of stationary solu- 
tions. For small pump values, the acoustic subharmonic 



2 



field is absent, and the homogeneous trivial solution is 
readily found as 



\X\=V, \Y\=0. 



(14) 



For higher values of the pump, the trivial solution be- 
comes unstable and the acoustic subharmonic field is 
switched-on. The corresponding values of the amplitudes 
are given by 



lATi-i, m , = l±VESpE£, (15) 

and the bifurcation occurs at a critical (threshold) pump 
given by 



Pth = 1, 



(16) 



which is independent of the value of the nonlinearity co- 
efficient rj. 

From (|15(l it also follows that the acoustic field shows 
bistability, i.e., both solutions coexist for pump values 
between Vb and Vth , where Vb corresponds to the turning 
point of the solution and is given by 



Hopf instability A = ±iu> is given by c^{c\C2 — C3) — c 2 C4 = 
0, or 

= 4af (7 - 1) \Yf + 2 V 2 (7 + if \Y\ 6 + (19) 
V 2 (4 7 3 - 3 7 - 3) |!f + 
2(7 + l) 2 (2 7 +l)|y| 2 -(2 7 +l) 2 . 

The frequency of the oscillations at the bifurcation 
point, given by the imaginary part of the eigenvalue, is 
found by substituting A = iui in the polynomial. In terms 
of the acoustic intensity given by the solutions of (|19fl the 
frequency reads 
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In Fig.l the domain of existence of Hopf bifurcations 
(the solutions of Eq. Q19J1) is represented in the plane < 
T), \ Y\ >, for different values of the relative loss parameter 
7 = 0.01,0.1,0.25 and 0.5. The stationary solutions are 
unstable at the left of each curve, i.e. for small values 
of i). From Fig. 1 also follows that, for a given value of 
the relative loss parameter 7 , a maximum value of r\ is 
required for the existence of dynamic solutions. 



1 +7? 2 



(17) 



We consider next the stability of the homogeneous so- 
lutions l|15l) . by means of the well known linear stabil- 
ity analysis technique. Substituting in Eqs. (|12|l and 
their complex conjugate a perturbed solution in the form 
Xi (t) =Xi + 5xi (t) , where Xi IS Si vector with the particu- 
lar stationary values, and linearizing the resulting equa- 
tions around the small perturbations, one obtains that 
Sxi ~ e At , where A are the eigenvalues of the stability 
matrix that relates the vector of the perturbations with 
their temporal derivatives. The instability of the solution 
is determined by the existence of positive real parts of the 
roots A of the fourth-order characteristic polynomial 
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where 



ci = 2(1 + 7), 

C2 = 4 7 + i + |y| 2 (77 2 |r| 2 -4 7 ) , 
C3 = 2 7 (i-2|r| 2 (1 + 7-^ |f| 2 )), 

C4 = 4 7 2 |11 2 (-l + |Y| 2 (l + 77 2 )), 

with |F| 2 given by Eq. ifHfl) . Since we are interested 
in the existence of dynamical behaviour, we search for a 
pair of complex conjugate eigenvalues, which denote the 
occurrence of a Hopf (oscillatory) instability. Following 
the Hurwitz criterion (see, e.g. 0), the condition for 



FIG. 1: Domain of existence of oscillatory solutions in the 
plane < 77, \Y\ >, for different loss parameters. Stationary 
solutions are unstable at the left of each curve. 



IV. SELF-PULSING DYNAMICS NEAR HOMO- 
CLINIC BIFURCATIONS 

In this section the numerical integration of Eqs. 1|12|) 
is performed in order to demonstrate the existence of dy- 
namical solutions. For typical experimental conditions, 
7q > 7(j, and consequently 7 < 1. We take, according to 



3 



[4|, the value 7 = 0.1 for numerical integration, although 
we note that similar results are obtained for a smaller loss 
ratio. For this value, the analysis of the previous section 
(see also Fig. 1, line c) shows that the solutions display 
temporal dynamics when the nonlinearity coefficient take 
values in the range < rj < 0.82. 
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FIG. 2: Bifurcation diagram of the ultrasonic field for 7 = 0.1 
and r) = 0.7. Open circles represent maximum and minimum 
amplitudes of oscillatory solutions. 



FIG. 3: (a) Time dependence of the intensity of the ultrasonic 
field for V = 1.35. The rest of parameters as in Fig. 2. (b) 
The corresponding phase portrait. 



in the form of periodic pulses or spikes, whose temporal 
separation (period) depends on the value of the pump 
amplitude, and increases with it. Close to the emis- 
sion threshold the time between pulses (interspike period) 
tends to infinity. An example of the temporal profile of 
the field amplitude in the spiking regime, and the corre- 
sponding phase portraits for the pump value V = 0.95 is 
shown in Fig. 4(a) and (b) respectively. 



In Fig. 2 the bifurcation diagram as computed numer- 
ically for 7 = 0.1 and r\ = 0.7 is shown. Dashed lines 
correspond to the analytical solutions given by Eq. (|15fl , 
and the existence of the backward (inverted) Hopf bifur- 
cation of the upper branch predicted by the theory of 
the previous section is observed at V ~ 1.87. Dynamical 
states exist for pump values below this critical point (see 
also Fig. 1). For higher pump values, the solution is sta- 
tionary, and the dynamics decay to a fixed point. The 
upper and lower branches with open circles in Fig. 2 cor- 
respond to maximum and minimum values of ultrasonic 
amplitudes in the oscillating regime. Numerics also show 
the switch-off of the ultrasonic field when pump is de- 
creased down to a given critical value. This means that, 
despite bistability, the trivial (zero) solution is a globally 
attracting solution, which have important consequences 
on the dynamics of the system. 

Similar diagrams are obtained for any value of ij in the 
range given above. For small values of the nonlinearity 
parameter ry, the solutions are quasi-harmonic in time re- 
gardless the pump value. However, for larger values of 
r] close to limiting value 0.82 required for the existence 
of dynamical states, the solutions show a qualitatively 
different temporal behaviour, not predicted by the linear 
stability analysis which in fact is valid only near the bi- 
furcation point. Next we report the numerical results for 
the particular value 77 = 0.7 corresponding to Fig. 2. 

In Fig. 3(a) the amplitude of the ultrasonic field is 
plotted as a funcion of time for a pump value oiV = 1.35 
far from the bifurcation point at Vth = 1- The corre- 
sponding phase portrait is shown in Fig. 3(b), where the 
existence of a stable limit cycle is observed. However, 
when decresing the pump, the ultrasonic field is emmited 
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FIG. 4: (a) Time dependence of the intensity of the ultrasonic 
field for V = 0.95. The rest of parameters as in Fig. 2. (b) 
The corresponding phase portrait. The homoclinic point is 
represented by the open circle. 

This behavior denotes the existence of a global bifu- 
raction of the solutions for the parameters considered. 
As the pump approaches the critical value V c = 0.904, 
the amplitude of the lower branch of the oscillating so- 
lutions (open circles in Fig. 2) approaches the lower 
branch of the stationary solutions (short-dashed line in 
Fig. 2), which as follows from the linear stability analy- 
sis corresponds to a saddle point. The limit cycle at this 
point connects with the saddle point, and denegerates 
in a saddle-loop or homoclinic bifurcation. For pump 
values below V c the trayectory decays to a fixed point 
corresponding to the trivial solution. 

The period of the limit cycle near a saddle- loop bifur- 
cation is governed by a characteristic scaling law. Lin- 
earization of the dynamics around the saddle leads to the 
following expression for the period [|| 
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where V c — V measures the distance to the homoclinic 
bifurcation (which is assumed small) and A is the eigen- 
value in the unstable direction of the saddle point. 

Figure 5 shows the numerically evaluated period of the 
self-pulsed solutions in the whole range where dynamical 
states exists. Note the divergence in the period at V c — 
0.904. In order to check the homoclinic character of the 
bifurcation, the inset shows the linear fit of the period 
against \n(V c — V) for pump values close to V c . The slope 
of the linear fit is found to be 51.4, in good agreement 
with the linear stability result 1/A = 54.2, demonstrating 
the existence of the saddle-loop bifurcation. 



300 



200 

T 









250 
200 














150 








\ 100 








-5.0 -4.0 -3.0 
MP-Phom) 




I , I 


I , I 





1.0 1.2 1.4 1.6 1.8 

P 



FIG. 5: Period of the oscillations as a function of the pump. 
The insert shows the linear fit. 

Finally, when decreasing the pump below the saddle- 



loop bifurcation point, the subharmonic field is switched- 
off. Trivial solution acts then as a globally attracting 
point. 



V. CONCLUSION 

A model of parametric generation of ultrasound in a 
ferromagnetic material has been considered, and its sta- 
bility analyzed, revealing the existence of a Hopf bifurca- 
tion leading to self-pulsing dynamics. For selected val- 
ues of the decay rates and the nonlinearity parameter, 
the system also show a saddle-loop bifurcation which re- 
sults in a spiking regime in the emitted ultrasound, where 
the frequency can take arbitrary low frequencies or fir- 
ing rates. It has been demonstrated 0] that dynamical 
systems where saddle-loop bifurcations and stable fixed 
points coexist have also a property called excitability, 
characterized by (a) perturbations of the rest state 
beyond a certain threshold induce a large response before 
coming back to the rest state, and (b) there exist a refrac- 
tory time during which no further excitation is possible. 
This properties, which are characteristic in several bio- 
logical problems (e.g. the behaviour of action potentials 
m neurons 

0) 

and some laser models |f|, are being 
investigated in the context of the acoustical system pro- 
posed in this paper. 



Acknowledgment 

The work was financially supported by the CICYT of 
the Spanish Government, under the project BFM2002- 
04369-C04-04. 



[1] A. Brysev, L. Krutyansky, P. Pernod and V. Preo- 

brazhensky, Appl. Phys. Lett. 76, 3133-3135 (2000). 
[2] L. Krutyansky, P. Pernod, A. Brysev, F. Bunkin and 

V. Preobrazhenshy, IEEE Trans. Ultrason., Ferroelect., 

Freq. Contr. 49, 409-414 (2002). 
[3] A. Brysev, L. Krutyansky and V. Preobrazhensky, Phys. 

Uspekhi 41, 793-805 (1998). 
[4] A. Brysev, P. Pernod and V. Preobrazhensky, Ultrasonics 

38, 834-837 (2000). 
[5] V.N. Streltsov, BRAS Physics Supplement 61, 228-230 

(1997). 

[6] V. Preobrazhensky, Jpn. J. App. Phys. 32, 2247-2251 



(1993). 

[7] A. Abd-Alla and G. Maugin, J. Acoust. Soc. Am. 82, 

1746-1752 (1987). 
[8] H. Haken, Synergetics, Springer- Verlag, Berlin (1983) 
[9] P. Gaspard, J. Chem. Phys. 94, 1 (1990) 
[10] E.M. Izhikevich, Int. J. Bif. Chaos 10, 1171 (2000) 
[11] J.D. Murray, Mathematical biology. Springer (1990) 
[12] A.L. Hodking and A.F. Huxley, J. Physiol. 17, 500 
(1952). 

[13] B. Krauskopf, K. Schneider, J. Sieber, S. Wieczorek and 
M. Wolfrum, Opt. Commum. 215, 367-379 (2003) 



5 



